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ABSTRACT 

Magnetohydrodynamic (MHD) simulations have been used to study disk accretion to a rotating mag- 
netized star with an aligned dipole moment. Quiescent initial conditions were developed in order to avoid 
j^j". the fast initial evolution seen in earlier studies. A set of simulations was performed for different stellar 

magnetic moments and rotation rates. Simulations have shown that the disk structure is significantly 
changed inside a radius r{, r where magnetic braking is significant. In this region the disk is strongly 
inhomogeneous. Radial accretion of matter slows as it approaches the area of strong magnetic field and 
a dense ring and funnel flow form at the magnetospheric radius r m where the magnetic pressure is equal 
, to the total, kinetic plus thermal, pressure of the matter. 

Funnel flows (FF), where the disk matter moves away from the disk plane and flows along the stellar 
magnetic field, are found to be stable features during many rotations of the disk. The dominant force 
driving matter into the FF is the pressure gradient force, while gravitational force accelerates it as it 
approaches the star. The magnetic force is much smaller than the other forces. The funnel flow is found 
to be strongly sub-Alfvenic everywhere. The FF is subsonic close to the disk, but it becomes supersonic 
■ well above the disk. Matter reaches the star with a velocity close to that of free-fall. 

' Angular momentum is transported to the star dominantly by the magnetic field. In the disk the 

• transport of angular momentum is mainly by the matter, but closer to the star the matter transfers its 

angular momentum to the magnetic field and the magnetic field is dominant in transporting angular 
momentum to the surface of the star. For slowly rotating stars we observed that magnetic braking leads 
to the deceleration of the inner regions of the disk and the star spins up. For a rapidly rotating star, the 
' inner regions of the disk rotate with a super-Keplerian velocity, and the star spins-down. The average 

torque is found to be zero when the corotation radius r cor w 1.5r m . 
^ ■ The evolution of the magnetic field in the corona of the disk depends on the ratio of magnetic to 

matter energies in the corona and in the disk. Most of the simulations were performed in the regime of 
a relatively dense corona where the matter energy density was larger than the magnetic energy density. 
' In this case the coronal magnetic field gradually opens but the velocity and density of outflowing matter 

are small. In a test case where a significant part of the corona was in the field dominated regime, more 
dramatic opening of the magnetic field was observed with the formation of magneto-centrifugally driven 
outflows. 

. i-h , Numerical applications of our simulation results are made to T Tauri stars. We conclude that our 

quasi-stationary simulations correspond to the classical T Tauri stage of evolution. Our results are also 
| relevant to cataclysmic variables and magnetized neutron stars in X-ray binaries. 

Subject headings: accretion, dipole — plasmas — magnetic fields — stars: magnetic fields — X-rays: 
stars 



1. INTRODUCTION 

Disk accretion to a rotating magnetized star is impor- 
tant in a number of astrophysical objects, including T 
Tauri stars (Edwards et al. 1994), cataclysmic variables 
(e.g., Warner 1995), and X-ray pulsars (e.g., Bildsten et 
al. 1997). 

The accretion of matter to a rotating star with a dipole 
magnetic field is a complicated problem still only partially 
solved. The important questions which need to be an- 
swered include: (1) What is the structure of the disk near 
the magnetized star? (2) Where is the inner radius of the 
disk? (3) What is the nature of the funnel flows (FF)? For 
example, which force is dominant in lifting matter to the 



funnel flow? (4) How is the accretion rate influenced by the 
star's magnetic moment /z and angular velocity 17*? (5) 
What is the mechanism of angular momentum transport 
between the star and the disk? What determines whether 
star spins-up or spins-down? (6) What are the necessary 
conditions for magneto-centrifugally driven outflows from 
the disk and/or the star? 

Many of these questions have been investigated analyt- 
ically, but the conclusions reached by different authors of- 
ten differ because the simplifying assumptions are differ- 
ent. For example, regarding question (2), some authors 
conclude that the disk should be disrupted in the region 
where magnetic and matter stresses are comparable (e.g., 
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Fig. 1. — Initial conditions of Types I (left panel) and II (right panel). The gray-scale and numbers show density 
distribution. 



Pringle and Rees 1972 - hereafter PR72; Davidson & Os- 
triker 1973; Lamb, Pethick & Pines 1973; Ghosh, Lamb, 
& Pethick 1977; Scharlemann 1978; Ghosh & Lamb 1979 
a,b - hereafter GL79a,b; Camcnzind 1990; Konigl 1991; 
Shu it et al. 1994 - hereafter S94). Others, argue that 
the inner radius of the disk should be farther away, at the 
corotation radius, because the inner regions of the disk 
are disrupted by magnetic braking (Ostriker & Shu 1995 
- hereafter OS95; Branderburg & Campbell 1998; Elstner 
& Riidiger 2000). 

Question (3) is investigated in only a few papers which 
use the Bernoulli integral (Lamb, Pethick & Pines 1977; 
Paatz & Camenzind 1996 - hereafter PC96; Li & Wil- 
son 1999; and Koldoba et al. 2002). The authors agree 
that the flow should become supersonic (and slow magne- 
tosonic) just above the disk. On the other hand opinions 
differ regarding the driving force which pushes matter up 
into the FF. Li and Wilson (1999) (see also Li et al. 1996) 
propose that the twisting of the magnetic field near the 
base of the FF should be very large, 7^ = \B^\lB p >> 1 
and the magnetic force should be the main one lifting mat- 
ter to the FF. Here, is the toroidal component of the 
magnetic field and B p is the poloidal component. Other 
groups (e.g., Lamb et al. 1977; PC96; Koldoba et al. 2002) 
argue that the FF should be super-Alvenic so that the 
twisting of the field is small so that the magnetic force is 
also small. In numerical simulations by MS97, H97 and 
GBW99 the magnetospheric accretion was reported, but 
no clear evidence of funnel flows was presented and no 
analysis of FFs was performed. A significant part of this 
paper is devoted to the FFs. Another important issue 
which has been discussed over the past 30 years is ques- 
tion (5) concerning the transport of angular momentum 
between the disk and the star. What determines the sign 
of the torque on the star? In early papers it was sup- 
posed that a star can only be spun-up because matter in a 
Keplerian disk brings positive angular momentum to the 
star (e.g., PR72). Later, it was recognized that a star can 
be spun-down due to the part of the star's magnetic flux 
which passes through the disk outside of the corotation ra- 
dius (GL79b, Wang 1995). Recently, the idea of "torque- 
less" accretion was proposed where mass but not angular 
momentum is transported to a star (e.g. S94; OS95; Li, 
Wickramasinghe & Riidiger 1996; Li & Wickramasinghe 
1997). Wang (1997) presented arguments against this idea, 
but this still remains an open question. 

Question (6) regarding magncto-ccntrifugally driven 
outflows from the disk has been discussed by a number of 
authors (e.g., Camenzind 1990; Konigl 1991; S94; OS95; 
Lovelace, Romanova & Bisnovatyi-Kogan 1995 - hereafter 
LRBK95; Fcndt, Camenzind & Appl 1995; PC96; Good- 
son & Winglee 1999; Bardou & Heyvaerts 1999; Agapi- 
tou & Papaloizou 2000). For example, S94 proposed that 
poloidal magnetic flux accumulates near the corotation 
radius and magnetic winds should blow from this point 
(X-point). LRBK95 proposed that wind may form from 
the entire region of the disk outside the corotation radius 
where the magnetic field threading the disk is open. 



Analytical investigations of disk accretion to a magne- 
tized star arc of course limited by the different assump- 
tions made. For this reason robust 2D and 3D simulations 
are essential to further the understanding of the different 
phenomena. By robust we mean that the result should 
not depend on initial conditions, boundary conditions, on 
grid resolution, and other artificial factors. Several 2D 
MHD simulation studies have been made with different 
initial conditions aimed at disk/star outflows. In an early 
work Hayashi, Shibata, and Matsumoto (1996) (hereafter 
HSM96) investigated the interaction of a non-rotating star 
with a Keplerian accretion disk and observed the opening 
of magnetic field lines which initially thread both the star 
and the disk. They found single event outflows and the 
corresponding inward collapse of the disk on a dynamical 
time-scale (less than one period of rotation of the inner ra- 
dius of the disk) with the radial velocity of the disk close 
to free-fall. This fast evolution was the result of the mag- 
netic braking of the disk by the magnetic field linking the 
disk and non-rotating star through a conducting corona. 
This explosive behavior may correspond to some episodic 
accretion events of actual systems. However, it is impor- 
tant to investigate the possible quiescent behavior of the 
disk-star systems. 

Miller and Stone (1997; hereafter MS97) investigated 
disk-star interaction for different geometries and stellar 
magnetic fields using the resistive ZEUS code. MS97 ro- 
tated the corona - which occupies all the space between 
the star and the disk - with the rotation rate of the star. 
This decreased the initial magnetic braking (compared to 
HSM96), so that they were able to perform simulations for 
several periods of rotation of the inner radius of the disk. 
In cases with a relatively weak magnetic field, they got re- 
sults similar to those of HSM96. They also found the disk 
collapsing to the star with velocity ~ 0.5w//, the open- 
ing of magnetic field lines, and outflows of matter from 
the disk. However, in the case of a strong magnetic field, 
particularly in the case which included a uniform homo- 
geneous vertical magnetic field threading the disk, they 
observed diminished outflows. Instead the matter flowed 
around the magnetosphere to the star. Similar results were 
obtained by Hirose et al. (1997; hereafter H97). 

Goodson, Winglee and Bohm (1997) (hereafter GWB97) 
and Goodson, B6hm and Winglee (1999) (hereafter 
GBW99) did much longer simulations in very large simu- 
lation regions. They observed quasi-periodic matter out- 
bursts associated with the quasi-periodic opening of mag- 
netic field lines and matter accretion to the star. The den- 
sity in the corona was chosen to decrease in a special way, 
1/i? 4 , so that the Alfven speed oc B/^/p oc 1/R decreases 
gradually. This is favorable for the opening of magnetic 
field lines, and for the generation and propagation of out- 
flows. GWB97 and GBW99 do not investigate cases where 
the density falls off more slowly with distance. 

Fendt and Elstner (1999, 2000 - hereafter FE99, FE00) 
investigated disk-star interaction for thousands of rota- 
tions of the inner radius of the disk, and observed the 
opening of magnetic field lines and outflows. However, 
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Fig. 2.— Evolution of a disk with no magnetic field for initial conditions of Type I (panels a, b, c, d) and Type II (panels 
e, f, g, h). Panels a & e show initial density distribution. Other panels show density distribution for different viscosities, 
a = 0, 0.01, and 0.02 at T — 50. The initial density for initial conditions of Type II is about 1.5 times smaller than for 
Type I. Consequently the density scale is slightly different in these cases. The legend for p is chosen to show the density 
distribution in the main part of the disk. It does not reflect the low density in the corona where p m in = 0.003 and the 
highest densities of matter in the disk near the star where p m ax = 18 at the panel c, p m ax — 12 at the panel d, p m ax — 2.1 
at the panel g, p m ax = 4.2 at the panel h. 



Fig. 3. — Evolution of the disk and the poloidal magnetic field in the largest region studied, R max = 55. The density 
(gray scale background) varies from a minimum value p — 0.003 in the corona to a maximum value p = 2.6 in the disk. 
The field lines are labeled by their magnetic flux \& values which change from 0.0009 to 4.9. 



Fig. 4. — Evolution of the disk and the magnetic field in a medium size region R m ax — 7. The configurations at 
T = 0, 10, 20, 30, 40, &50 are shown. The density (gray scale background) changes from maximum value p = 2.6 in the 
disk to the minimum value p = 0.003 in the corona. Contour levels of "J, which label poloidal field lines, change from 4.9 
to 0.006. The field lines of the strongest magnetic field near the star are not included in order to make the inner structure 
of the disk visible. 



they treated the disk as a boundary condition so that they 
could not take into account the back reaction of the disk 
on the stellar magnetic field. Furthermore, the actual out- 
flow of matter from the disk to the corona may be different 
from that assumed. 

The above - mentioned simulation studies show that 
outflows appear either in very non-stationary situations 
(HSM96, MS97) or for a special distribution of coro- 
nal density and very fast rotation of the star (GWB97, 
GBW99). None of the papers give answer to questions 
(l)-(5). Also, it is not clear whether or not outflows exist 
for more quiescent initial conditions, and for cases where 
the coronal density falls off slowly with distance. 

In this paper we investigate disk accretion to a rotat- 
ing magnetized star and the associated funnel flows. We 
start from initial conditions which give us the possibility to 
significantly reduce the initial magnetic braking between 
the disk and the corona. This allows us to investigate the 
disk-star interaction and funnel flows for long times and 
to consider questions (l)-(5) in detail. 

In §2 we describe the numerical model, including initial 
and boundary conditions. We also discuss the evolution 
of the disk without a magnetic field. In §3 we describe in 
detail the disk-star interaction for the case of slowly ro- 
tating stars and in §4 cases of fast rotating stars. In §5 
we consider the dependence of disk-star interaction on the 
magnetic moment p. In §6 we analyze the physics of FFs. 
In §7 we consider the possibility of outflows. In §8 we ap- 
ply our simulation results to T Tauri stars. In §9 give the 
conclusions from this work. 

2. NUMERICAL SIMULATIONS 

To investigate the disk-star interaction, we numerically 
solve the MHD equations, 



dp 
dt 
dB 



+ V-(pv) = (1), 



dt 



-V x (v x B) = (3), 



dt 

d(pS) 



v • r = p g 



(2) 



+V • (pvS) - . (4) 



dt 

Here, S is entropy, T jk = p5 jk + pvjv k -r ^ u Jk 
BjB k )/(4:Tr) is the stress tensor, g = — V$ is the gravita- 
tional acceleration, and $ = —GM/R is the gravitational 



(B 2 S jk /2- 



potential of the central object. We take the equation of 
state to be S = p/p 1 , where 7 = 5/3 for most of this 
work. We solve equations (l)-(4) in a spherical coordinate 
system (R, 9, <p). 

In order to model the slow accretion of matter in the disk 
we modified equations (2) and (4) by including a small vis- 
cous stress Gjk in Tjk- Because the flow is dominantly in 
the azimuthal direction the main viscous stress acts in this 
direction. Thus the important viscous stress components 
are the (R<t>) and (0(f)) terms, ctr0 = —vpRsinOdui/dR, 
aecf, = —vpsm9du)/d6, where u = v^/Rsind is angular 
velocity of the matter, and v is the kinematic viscosity. 
We adopt the a— model for the viscosity (Shakura & Sun- 
yaev 1973), where v = olc 2 s /VLk, with c s = (p/p) 1 / 2 the 
isothermal sound speed velocity and a < 1 is a dimen- 
sionless parameter, which was chosen to be a = 0.01 or 
a = 0.02 in most of simulations which gave a reasonable 
speed of matter flow to the central regions. We checked 
that for the grid sizes used the numerical viscosity is much 
lower than the considered a viscosity (see §2.5). 

Equations (1) - (4) were solved with a Godunov type 
numerical code developed and tested by Koldoba et al. 
(1992), Koldoba and Ustyugova (1994) and Ustyugova et 
al. (1995). This type of code has also been described 
by Ryu, Jones, and Frank (1995). A modified version of 
the code incorporating viscosity was tested in a number of 
hydrodynamic simulations for different values of a, from 
a = 0.01 to a = 0.1. For a ~ 0.01 - 0.05 we observed 
a steady flow of matter to the star with small velocities 
(see §2.5) corresponding approximately to the Shakura- 
Sunyaev accretion model. 

Our Godunov type code does not explicitly include 
magnetic diffusivity (unlike the codes of MS97, H97, & 
GBW99). However, matter diffuses across the magnetic 
field owing to the small numerical diffusivity v m which is 
of the order of the numerical viscosity v m <~ ARc s , where 
AR is the grid spacing. The estimated magnetic Reynolds 
number associated with the radial flow is Re m ~ 8Rv r /v m , 
where 5R is the characteristic scale. In the region 1 < 
R < 5, Re m ~ 3 owing to the small radial velocity v r . 
In actual accretion disks, both angular momentum trans- 
port and diffusion of the magnetic field, are probably due 



4 



Fig. 5. — Evolution of the density in the disk (color background) and poloidal magnetic field lines ('5 =const lines) (white 
lines) in case of Type I initial conditions in the region r < 6. The black solid line corresponds to (3 = 1. The density 
changes from 0.003 in the corona to 2.4 in the disk. Contours of \P are exponentially spaced between 0.2 and 1.2. 

Fig. 6.— Same as on Figure 4 but for initial conditions of Type II. At T = the inner radius of the disk is at r — r cor = 3. 
The ^ =const lines are exponentially spaced between 0.1 and 0.7. 



to turbulence (e.g., Bisnovatyi-Kogan & Ruzmaikin 1976; 
Balbus & Hawley 1991) and have similar order of magni- 
tude transport coefficients. In our simulations we observed 
that both the viscosity and the magnetic diffusivity are 
small and are of the same order of magnitide. Note, that 
in the <fi direction, the velocity of the flow is high, and the 
Reynolds number is also high, (i?e m )^ ~ 100, so that mag- 
netic field lines rotate with the disk in <p direction, and at 
the same time matter slowly diffuses in r— direction. This 
is important for understanding of evolution of magnetic 
field lines observed in the following sections. 

2.1. Reference Values and Dimensionless Units 

A reference value of distance is denoted Ro- The refer- 
ence value for velocity is taken to be v = (GM/i? ) 1 ^ 2 
which is the Keplerian velocity at Ro- The reference an- 
gular rotation rate is luq — vq/Rq, and the corresponding 
time-scale is to — Ro/ v o- For discussing our results we also 
use the rotation period at r = Ro, Po = 2nt - For a given 
magnetic field strength at Ro we can define a reference den- 
sity as po = Bq I Vq and reference pressure po = pqVq ■ A 
reference magnetic moment of the star is then po = BqRq. 
Thus, the calculated variables are: R 1 — R/Ro, t — t/to, 
p' = p/Po, v' = v/vq, B' = B/Bo, p' = p/po- We also 
use the dimensionless time, T — t/Po- Later, we omit 
the primes but note that any dimensionless variable can 
readily be converted to its physical value. 

Here, we give parameters for a typical T Tauri star. We 
take the mass and radius of the star to be M* = 0.8M Q and 
i?* = 2Rq. We use as the length scale Ro the initial inner 
radius of the disk, so that the dimensionless inner radius of 
the disk is Rd = 1 in most of simulations. We place the in- 
ner boundary (a star) at R m in = 0.35. Then, our reference 
length is R w i?*/0.35 w 4.0 x 10 11 cm, while the simu- 
lation region corresponds to R m ax = 55 x i?o ~ 1-47 a.u.. 

The reference velocity is v ~ 1.63 x 10 7 cm/s, and the 
corresponding time-scale is to ~ Ro/ v o ~ 2.45 x 10 4 s. 
The period of rotation of the inner radius of the disk is 
P ~ 1-78 days. 

Consider a magnetic field strength B* = 10 3 G at the 
star's surface. Then at R = Ro, the reference magnetic 
field is B = B*(R*/R ) 3 ~ 42.9 G and the reference 
magnetic moment is po = BqRo ~ 2.7 x 10 36 Gem 3 . 
The reference density is po — 6.9 x 10~ 12 g/cm 3 , or 
n = 6.93 x 10 12 cm~ 3 which is typical for T Tauri star 
disks. The reference mass accretion rate is Mo — Po v oRq ~ 
1.8 x 10 18 g/s w 2.8 x 10~ 7 M Q /year. The reference an- 
gular momentum flux is: L = po^o^o- 

2.2. The Grid 

We place the inner boundary of the simulation region - 
the "star" - at R m i n — 0.35. The size of external boundary 
depends on the grid (see below). The spherical grid was 



inhomogeneous in the R— direction. The inhomogeneity 
was such that cells at any distance R were approximately 
square with AR ~ RAO for fixed AO. This grid gives a 
high space resolution close to the star which is important 
in this problem. At the same time simulations in very 
large regions may be performed. The inhomogeneous grid 
is a smooth analog of the nested grids used by GBW97 
and GBW99. 

The angular grid was taken to have Ng = 51 in most 
cases, and Ng = 71 in cases with very strong magnetic 
field (sec §7). The number of points in the R— direction 
determined the size of the simulation region. The main 
runs were done with a large radial grid Nr = 150 which 
corresponds to R ma x = 55 ss 157R m i n . Simulations with 
a smaller grid Nr = 100 corresponding to R max — 10 ss 
35R m i n were done in a number of cases. The large size of 
the region was chosen in order to minimize the influence 
of external boundaries on processes in the vicinity of the 
star. We checked, however, that results of simulations in 
the larger region almost coincide with those in the smaller 
region. For the typical grid with Ng = 51, the smallest 
cell near the star corresponds to AR = 0.01 s=s 0.03i? m j n 
which is a resolution about 3.2 times higher than one in 
GBW99. The largest grid is at the external boundary, for 
Nr — 150 and is AR = 5AR min . These numbers show 
that a spherical stretched grid is advantageous for simula- 
tions of accretion flow to a star with a dipole field. 

2.3. Boundary Conditions 

At the inner boundary R = R m in, we apply "free" 
boundary conditions for the density, pressure, and en- 
tropy, dp/dR = 0, dp/ OR = 0, and dS/dR = 0. The 
inner boundary is treated as a rotating perfect conductor 
J7 = fiz, which is analogous to the case of a rotating disk 
discussed earlier by Ustyugova et al. (1995, 1999). In the 
reference frame rotating with the star, the flow velocity v 
is parallel to B at R = R m in ', that is, B x (v — fi x R) = 
at Rmin- We consider the "free" boundary condition for 
this velocity, <9(v — (I x R)n/dR = 0. These two con- 
ditions determine the direction and value of the velocity 
vector v, that is all three components ("Wr, t>e, v^). The 
boundary condition at R m in on the magnetic field has 
d{RB<j,)/dR = 0, while the poloidal components Br and 
Bg are derived from the magnetic flux function ty(R,0), 
where *f? at the boundary is derived from the frozen-in 
condition d^/dt + v p • W = 0. 

At the outer boundary R = R ma x, we took fixed bound- 
ary conditions for all variables for those cases when the 
simulation region was very large: R ma x — 55. In case of 
smaller region R max = 10 we took free boundary condi- 
tions on all the hydrodynamic variables and the B^ com- 
ponent of the magnetic field, if matter outflows from the 
region. For the magnetic field we took dB^/dR = and 
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d'i'/dt + v p • V\& = 0. If matter inflows to the region 
(vn(R max ) < 0), then we set this velocity equal to zero. 
Results are very similar at these boundary conditions. We 
assume reflection symmetry about the equatorial plane. 

2.4. Initial Conditions 

Here we present our new method for establishing quasi- 
cquilibrium initial conditions. The star has a magnetic 
dipole field with the magnetic moment fi aligned with 
the rotation axis z, with the star's rotation rate. The 
magnetic field is B = 3(/x • R)R/i? 5 — /i/i? 3 , with com- 
ponents in spherical coordinates Br = 2/xcos# / R 3 and 
Bq = psm9/R 3 . As initial conditions we set up a low- 
temperature Td disk with a high-temperature T c >• Td 
corona filling the remainder of the simulation region. The 
disk rotates with angular velocity close to the Keplerian 
value lu w Qk- The disk extends inward to the radius 
Rd where the ram pressure of the disk is equals to the 
magnetic pressure of the dipole, p + pv 2 = B 2 /8ir. 

In order to have approximately equilibrium initial con- 
ditions it is necessary that the corona initially be rotating. 
If the corona initially does not rotate (e.g., HSM96), or 
rotates with the angular velocity of the star (e.g. MS97), 
while the disk rotates with Keplerian velocity, then the 
coronal magnetic field lines threading the disk and the 
star lead to magnetic braking of the disk and subsequent 
fast accretion of the disk with a velocity close to free-fall. 
Furthermore, a strong discontinuity develops in the angu- 
lar velocity between the disk and the corona (except at 
the corotation radius). This leads to the generation of a 
strong Bfj, component, while the dragging of the magnetic 
field with the radial inflow of the disk leads to the gen- 
eration of a B r component. This twisting of the field is 
responsible for transient magneto-centrifugally driven out- 
flows (e.g HSM96; MS97; Hirose et al. 1997). Such fast 
evolution may correspond to periods of violent accretion in 
some strongly non-stationary systems. However, the star- 
disk systems may exist in more quiescent configurations 
which can be studied by using appropriate initial condi- 
tions. 

In order to avoid the strong discontinuity of lo and B$ on 
the boundary between the disk and the corona, we rotate 
the initial corona so that its angular velocity is a constant 
on a given cylindrical radius r = Rsin9 and equal to the 
Keplerian rotation rate of the disk. Thus the corona ro- 
tates with different angular velocities on cylinders with 
different radii r. For such initial conditions we observe 
much slower accretion in the disk with no dramatic out- 
flows. This distribution of u> in the corona also leads to the 
twisting of the magnetic field lines. Any coronal magnetic 
field line will be twisted along its length by the differential 
rotation of nearby layers of corona, but this twist is dis- 
tributed along the length of the field line. At the boundary 
between the disk and corona the twist is small. It is much 
smaller than the twist at the boundary between a rotat- 
ing disk and non-rotating corona. Differentially rotating 
corona also leads to the magnetic braking of the disk, but 
it is more gradual and does not lead to a rapid infall of the 



matter in the disk. Earlier, a differentially rotating corona 
was used as an initial condition in a study of the opening 
of magnetic loops threading a Keplerian disk, but the disk 
was considered as a boundary condition (Romanova et al. 
1998). 

Initially, we assume that B^ — everywhere. Also, we 
assume that the initial poloidal field is that of a dipole and 
that the plasma is force-free. Thus, we search for equilib- 
rium initial conditions for the disk-corona system taking 
into account gravity, pressure, and rotation. To find the 
initial equilibrium, we suppose that initially (t — 0) the 
matter is barotropic p = p(p). In this case one have (I) 
Vp/ p = VF, where F = J dp/ p; (2) the angular velocity w 
depends only on r = Rs'm9; (3) the centrifugal force has 
the potential $ c = — fZ [u(r')] 2 r'dr'; (4) the stationary 
Euler equation (v • V)v + p/p + V<f> = has the integral: 

F + <£> c + <S> = E = const, (5) 

where $ = —GM/R is the gravitational potential. 

We suppose that the initial angular velocity of the mat- 
ter in the equatorial plane is 

w(r,0-7r/2) - | {kGM/r3)K r > R ^ (6) 

where r = R sin 9 is the cylindrical radius, and k — const ~ 
I is included to take into account that the disk is slightly 
non-Keplerian (in our simulations k — I. Of). Thus, the 
centrifugal potential in the whole region is 

W " J "| kGM/r , r>R d . 

At the boundary between the disk and corona the initial 
pressure is taken to be po- The dependence of the density 
on pressure in the corona and in the disk can be expressed 
as 

o(o) - { ™ p/Tc ' P<Po (8) 
W> - \ mp/Td , p> Po , (S > 

where m is the mean particle mass approximately equal to 
m#/2. We consider that T d <C T c . In the corona, at low 
values of pressure and density (p < po,p < fhp /T c ) the 
plasma is hot T = T c , whereas in the disk, at high pres- 
sure and density (p > po,p > fhpo/Td) the plasma is at a 
much lower temperature, T — Td- Initially, the boundary 
between the disk and corona corresponds to the contact 
discontinuity, where p = po, p c = mp /T c in the corona, 
and pd = fhpo/Td in the disk. In this case 

p _ / (Tc/m)\n(p/p ) , p < po (Q , 
\ (Td/m)\n(p/p ) , p > p . {y > 

The initial disk-corona boundary in the equatorial plane 
is assumed to be at r = Rd- We assume values of po, T c , 
and Td- At the point (r = Rd, z = 0), p = po, F(p ) = 0. 
Thus, we find E = -GM/R d +kGM/R d = (k-l)GM/R d . 
Then, we can calculate potential + at any point (i?, 9) 



Fig. 7. — Radial distribution of angular velocity along the equatorial region of the disk for r < 2 at different times T. 
The solid line shows the Keplerian velocity. 
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Fig. 8.— Distribution of angular velocity in the disk and corona at different times T for the same case as Figures 4 and 
5. Only the inner region of the disk r < 4 is shown, u) changes from ui = 0.1 (white) to w = 1 (black). The solid lines are 
poloidal field lines (<J/ =const lines). 



and find F = E — ($ + <1> C ). Now, we can derive pressure 
p at any point of the region, 

J po exp (Fm/T c ) , p < Po nn x 
p -\ Po exp(Ffh/T d ) , p> Po ^ W > 

and find the distribution of density taking into account 
eq. (8). We considered two main types of initial condi- 
tions which are described below. 

2.4.1. Initial Conditions of Type I 

In the first type of initial conditions, we place the in- 
ner radius of the disk at r = Rd = 1. We determine the 
characteristic density and temperature in the corona to be 
p c = 0.01, T c — 1, and density in the disk pd = 1. Then we 
determine the pressure at the boundary between the disk 
and corona, po = p c T c = 0.01 and then derive tempera- 
ture in the disk, Td = po/pd = 0.01, so that the pressure is 
continuous across the boundary. Then we determine the 
magnetic field of the star so as to have the magnetic pres- 
sure equal to the stagnation pressure of the disk matter 
at R = 1. This gives p 2 /8nR^ w PdV 2 K which is solved 
for p. As matter in the disk flows inward, it forms a fun- 
nel flow at R i=a 1. We determine u(r) from equation (6) 
taking into account that Rd = 1. For r < Rd matter ro- 
tates with an angular velocity corresponding to %(r = 1), 
while at r > Rd the rotation is Keplerian. The surfaces 
of constant u> are cylinders. For such a distribution, there 
are no strong gradients of uu in the region. Then we de- 
rive the pressure and density distributions using equations 
(5)-(10). 

The left panel of Figure 1 shows the density distribu- 
tion in the disk and corona. For this example we took a 
relatively small region R m ax = 7. One can see that the 
density is high in the disk, and is 100 — 200 times smaller 
in the corona. For these initial conditions the FF starts 
close to the initial inner radius of the disk. Thus the disk 
provides a "reservoir" of matter for the FF. On the other 
hand there is the following disadvantage. In most cases, 
we rotate the star with an angular velocity f2* smaller than 
the initial angular velocity of matter above the star. For- 
tunately, the star is observed to control the rotation of the 
nearby regions of corona in a few time-steps. We also used 
initial conditions of Type II, as described below. 

2.4.2. Initial Conditions of Type II 

In a second type of initial conditions we place the in- 
ner radius of the disk R d at the corotation radius, r cor = 
(GM/nl) 1 / 3 . Thus for r < R d both the star and the 
corona rotate with the angular velocity Qk(t = Rd), while 



at larger distances the rotation is Keplerian. For this ini- 
tial condition, u> varies smoothly including the variation 
from the star and to the corona above it. 

An important aspect of the Type II initial conditions 
is that the initial inner radius of the disk is far from the 
magnctospherc. This allows us to check whether or not 
matter moves inward from the corotation radius. 

The left panel of Figure 1 shows the initial density dis- 
tribution for Type II initial conditions for r cor = 3. 

2.5. Hydrodynamic Evolution of the Disk and Corona 

First, we tested our initial conditions with no magnetic 
field. We observed that for both types of initial condi- 
tions the disk exists for more than 300 rotational periods 
Po- No high velocities appear in the disk or corona. Typ- 
ical velocities in the disk are Vd ~ 0.01 — 0001. However, 
our typical time of simulation of magnetic cases is shorter: 
50Po, so that we show the result of the evolution of the 
disk without a magnetic field at t < 50Pq- 

We considered three cases, with no viscosity, a = 0, 
and with small viscosities a = 0.01 and a = 0.02. Figure 
2 (a,b) shows that at a = and boundary conditions of 
Type I, the disk's inner edge stays at r = 1 during the 
whole period of evolution. The density in the rest of the 
disk slowly redistributes itself in such way that the external 
layers move gradually outward. The outer regions of the 
disk had excess angular momentum compared to a station- 
ary a-disk (Shakura & Sunyaev 1973), and consequently 
some matter moved outward carrying away extra angular 
momentum. When we included a small viscosity a = 0.01 
(see panel c of Figure 2), the outer regions of the disk also 
moved outward, but inside the shown region r < 10 matter 
moves inward with velocity v < 0.0001 — 0.001 for r > 6, 
and with larger velocity v ~ 0.003 — 0.006 for r < 6. 

For even higher viscosity a = 0.02 (panel d of Figure 2), 
the velocity of the inflow is v ~ 0.003 — 0.006 in the entire 
disk, r < 10. In both initial conditions we find that the 
a— viscosity gives a slow accretion and thus a good start- 
ing point for our MHD simulations. The typical accretion 
rate isM« 0.05 for a = 0.01 and M w 0.1 for a = 0.02. 

For the Type II initial conditions we observed that for 
a = the inner radius of the disk stays at r = 3 for T = 50. 
This shows that the numerical viscosity is very small (see 
panels e and / of Figure 2) compared with ana= 0.01 
viscosity. External regions moved outward with velocity 
v < 0.001, redistributing density and angular momentum. 
For a = 0.01 and 0.02, the disk behaved similarly to one 
for initial conditions of Type I: the equatorial regions ac- 
creted slowly inward with velocity v <~ 0.003 — 0.007 (pan- 
els g and h of Figure 2). The accretion rate is M w 0.02 
for a = 0.01 and M w 0.05 at a = 0.02. These tests 



Fig. 9. — The figure shows the time dependence of the mass accretion rate M and the angular momentum fluxes carried 
by matter L m and magnetic field L/ to the star for the case of a slowly rotating star fi* = 0.19 which corresponds to 
P* = 9.4 days). The corotation radius is r cor = 3. 
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Fig. 10. — Same as on Figure 9, but for initial conditions of Type II. 

Fig. 11. — Results of simulations of disk-star interaction for the case of a rapidly rotating star, fi* = 1 and r cor = 1, 
which corresponds to P* w 1.8 days for the T Tauri parameters of §2.2. The density (background) and poloidal magnetic 
field lines or * =const lines (solid lines) are shown. The \E" contours are exponentially spaced between 0.15 and 1. 



Fig. 12.— Same as on Figure 9, but f2» = 1 and r cor = 1 (P* 

indicate that our hydrodynamic disk and corona are ap- 
propriate for subsequent MHD simulations of the disk-star 
interaction. 

Note, that subsequently we observed that accretion rate 
to the star in the magnetic case is larger than in hydro case 
(see Figures 9 and 10, where accretion rate is M w 0.2 — 0.4 
at initial conditions of Type I, and M w 0.5 — 1.3 at ini- 
tial conditions of Type II, which is connected with small 
magnetic braking of matter in the disk, because initial con- 
ditions give only approximate but not exact initial equi- 
librium. 

3. DISK-STAR INTERACTION FOR A SLOWLY ROTATING 

STAR 

Here, we discuss simulations of accretion to a slowly ro- 
tating star with fi* = 0.19 (r cor = 3), which corresponds 
to a period P* = Po/O* « 9.4 days for a typical T Tauri 
star (see §2.2). The magnetic moment is p = 1.06, which 
corresponds to a magnetic field at the surface of the star 
P* w 1100 G. We used initial conditions of Types I and 
II to check for the possible dependence of the evolution on 
initial conditions. The simulations were performed with 
the grid Nr x Ne = 150 x 51 which has R ma x = 55. 

3.1. Simulations with Type I Initial Conditions: 
Structure of the Disk 

First, we describe results of simulations with Type I ini- 
tial conditions. We observed that matter started to move 
slowly both in the disk and corona with velocities much 
smaller than free-fall. No dramatic collapse of the disk or 
fast outflows were observed owing to our quasi-equilibrium 
initial conditions. 

Figure 3 shows the entire simulation region at different 
times. The outer part of the disk changes very little during 
the simulation time T < 50. The magnetic field lines twist 
and become partially open. Viewed on this large scale, the 
disk is relatively flared, but only the much smaller region 
of the disk R < 10 is important for disk-star interactions. 

Figure 4 shows the evolution in the smaller region R < 7. 
One can see, that the inner region of the disk is influenced 
by the magnetic field of the star. For T s=s 50 the accre- 
tion rate due to the a— viscosity increased and magnetic 
flux was accumulated in the central region. We stopped 
the simulations at this point, because later the accretion 
rate increased, the inner radius of the disk moved closer 
to the star, and the case became less interesting for the 
investigation of the disk-star interaction. 

Figure 5 shows the evolution of the density and poloidal 
magnetic field in the region r < 6. The interaction of the 
disk with the magnetosphere of the star led to the recon- 
struction of the disk inside the radius which we call the 
"braking" radius ?v. In this case 7V w 3. This radius is 



= 1.8 days for the T Tauri parameters of §2.2). 

an approximate analog of the radius of the outer transi- 
tion zone in the theory of GL79b, which divides areas of 
an undisturbed viscous disk and inner regions, influenced 
by a magnetic field. In our case, the disk for r > r\, r is also 
disturbed by the star's field, but matter still moves to this 
region due to viscosity. At smaller radii r <; 7V, matter is 
magnetically braked. By coincidence rt r w r cor . The in- 
ner regions of the disk become magnetically braked after 
several rotation periods p). In the region of the magnetic 
braking the density in the disk is 2— 4 times smaller than in 
the disk without a magnetic field. The inflow speed of the 
matter increases to v <~ 0.01 — 0.06. and matter constantly 
moves in the direction of the star. The flow towards the 
star stops when the matter stress becomes comparable to 
the magnetic pressure of the dipole, B 2 /8n = p + pv 2 . 
Note that the main term in the matter stress, pv 2 ^, is due 
to the disk's rotation. We introduce the plasma parame- 
ter, (3 = (p + pv 2 )/{B 2 /8-k), which indicates the relative 
importance of the matter and magnetic stresses. Also, we 
introduce the magnctosphcric radius, 



m [8w(p + pv 2 )}V e ' ^ 

which is the radius at z — where (3 = 1. 

Figure 5 shows that matter accumulates in the region 
r ~ r m forming a dense ring. The density in the ring is 
^2 — 3 times larger than the initial density at the inner 
radius of the disk. At T = 50, when more matter comes 
from the disk, the density in the ring becomes ~ 10 times 
larger. Accumulation of matter near r m and formation of 
a ring in the inner regions of the disk was also reported by 
GWB97, GBW99, and MS97, but these simulations were 
for non-stationary conditions. 

For r <J r m matter moves upward, away from the equa- 
torial plane and goes into the funnel flow. The FF is 
formed after about one period of rotation (T <~ 1) and 
is a stable feature during the length of the run, T = 50. 
Figure 5 shows the (3=1 line (black line) . The base of the 
FF coincides with r = r m ; that is FF starts in the area 
where the magnetic stress is equal to the matter stress. 
This is observed in many other cases with different pa- 
rameters (see §4 & §5). This is in accord with theories 
developed in 1970's (PR72, GL79a,b). We investigate the 
physics of FFs in §6. 

Matter from the outer region of the disk r > rv tends 
to accrete inward. Often accretion rate is enhanced in the 
region where magnetic field has larger radial component 
B r and hence stronger magnetic braking. 

The accretion rate of matter from the region of the disk 
outside the ring varies with time. Consequently, the den- 
sity in the inner ring and in the funnel flows also varies 
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with time (Figure 5). When more matter accumulates in 
the ring, then the base of the FF is closer to the star, 
and the magnetospheric radius is smaller. For example, at 
T = 10 (panel b of Figure 5), r m = 0.92; at T = 30 (panel 
c) r m = 1.2. Subsequently, for T > 50, even more matter 
came from the disk owing to viscosity, and the FF moved 
even closer to the star. A stationary state for T > 50 may 
exist but at smaller values of a, a < 0.02. 

3.2. Simulations with Type II Initial Conditions 

For Type II initial conditions, the inner radius of the 
disk is taken to be the corotation radius, r cor = 3. The 
viscosity parameter is a = 0.02, the same as in the above- 
mentioned simulations. We were able of observe the evo- 
lution for more than 80 rotations of the inner radius of the 
disk. Figure 6 shows the evolution of the inner regions of 
the disk. We observed that the funnel flow did not form at 
the corotation radius as predicted in some theories (e.g., 
OS95). Instead, the disk slowly moves inward towards the 
star owing to the viscosity (see Figure 6). When the disk 
reached conditions where /3»1, the accreting matter goes 
up into the funnel flow. The FF formed after T ss 10 
rotations. At longer times, the structure of the disk and 
the FF are similar to those found for Type I initial con- 
ditions (Figure 5). But one difference is that for Type 
II initial conditions the inward motion of the inner edge 
of the disk drags the bases of magnetic loops inward, so 
that loops of the dipole are inclined to the disk. Such a 
situation is possible when the viscosity is larger than the 
diffusivity and magnetic flux in the disk accumulates closer 
to the star (e.g., S94, OS95), but the rest of the loop in 
the corona doesn't move so fast, so that the loops are in- 
clined to the disk, and the B r component is significant. 
When the magnetic field is inclined away from the axis 
of the disk, the magnetic force is much stronger than in 
case with no inclination (Lovelace, Berk & Contopoulos 
1991, Oyed & Pudritz 1997; FE99, FE00). We observed 
opening of the magnetic field lines close to the magneto- 
sphere. The opening was followed by reconnection, and 
subsequent opening. The process was quasi-periodic with 
a quasi-period T w 3.5 days. We also observed quasi- 
periodic oscillations of the inner radius of the disk between 
radii: 0.7 < r m < 1.2. The nature of the oscillations 
is the following: (1) First, the loop is stretched and the 
magnetic field "blocks" the path of matter to the FF. (2) 
Next, reconnection releases the stress, and matter accretes 
through the FF. (3) The new loop starts to stretch, and 
the process repeats. Figure 10 shows the quasi-periodicity 
of the accretion rate. These oscillations resemble those 
observed by GWB97 and GBW99, but on a much smaller 
scale. H97 and MS97 also observed enhanced accretion af- 
ter the reconnection of elongated magnetic field lines near 
the magnetosphere. This work supports the idea of Aly 
& Kuijpers (1990) who proposed this mechanism of quasi- 
periodic oscillations. 

We tested the dependence of the results on the viscos- 
ity for both types of initial conditions for a = 0.01 and 
a = 0. We obtained similar results, but the FF formed at 
slightly larger distances r as a result of smaller accretion 



rate. Simulations with initial conditions of Types I and II 
have shown examples of possible situations which can be 
realized in nature. In the remainder of the paper we use 
the Type I initial conditions, because the FF starts earlier 
in this case. 

3.3. Angular Momentum Transport 

How is angular momentum transported between the disk 
and the star? Does the star spin-up or spin-down? In this 
section we analyze these questions using the run described 
in §3.1. 

3.3.1. Angular Velocity of the Disk 

We observed that the magnetic field lines which are re- 
sponsible for the partial destruction of the disk at radii 
r < rv are also involved in the transport of angular mo- 
mentum between the star and the disk. The magnetic filed 
lines connecting the star and the disk decrease the disk's 
angular velocity in a significant part of the braking region. 
Figure 7 shows the radial distribution of the angular ve- 
locities at different times. One can see that the angular 
velocity of the disk is significantly smaller than the Keple- 
rian value in a broad region r < 1.3 — 1.8. We introduce the 
radius r*, such that for r < r^ the magnetic force on the 
disk is strong enough to cause a significant deviation from 
Keplerian rotation. The radius depends strongly on 
the accretion rate and varies between r$ w 1.8 at T = 30, 
(when the inner disk is of low density) and r^ w 1 at 
T = 10 (when the inner disk is dense). 

Figure 8 shows the spatial distribution of the angular 
velocity, oj(R,8). One can see that for r < r\, r rts 3 the 
regions of constant angular velocity uj(R, 0) = const, al- 
most coincide with the regions of constant magnetic flux 
Thus, the magnetic field line crossing the disk at radius 
r rotates in the corona with angular velocity u>(r). This 
is because most of the corona is matter-dominated, f3 > 1 
(see black line in Figure 5), and the disk determines the 
rotation of the magnetic field lines. Magnetic field lines 
which cross the disk in the magnetically-dominated area 
(in the FF) rotate with the angular velocity of the star. 

3.3.2. Spin Evolution of the Star 

An important question is the angular momentum flux to 
the star. To answer this question, we calculated the flux 
of angular momentum (about the z— axis) carried by the 
matter L m and that carried by the 'twist' of the magnetic 
field Lf, 

L m = J dS • prv^Vp , L f = - J dS • rB^Bp/An , 

evaluated at the surface of the star. In addition, we cal- 
culated accretion rate to the star, M — J dS • pv. 

Figure 9 shows the different fluxes as a function of time. 
Most of the angular momentum flux at the star is carried 
by the magnetic field; matter carries only <~ 1% of the to- 
tal flux. The flux Lf is positive so that it acts to spin-up 



Fig. 13.— Same as on Figure 9, but fi* = 0.54 and r cor — 1.5 (P* w 3.3 days for the T Tauri parameters of §2.2). 
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Fig. 14. — Disk-star interaction in case of "torqueless" accretion at fi* = 0.54 and r cor = 1.5 (P* = 3.3 days for the T 
Tauri parameters of §2.2). Contour lines of magnetic flux ^ are exponentially spaced between 0.07 and 0.7. These values 
of ^ are smaller than those in Figure 11 because density at the inner edge of the disk is smaller and disk stops at smaller 
values of 4'. 



the star. The flux Lf correlates with M because incom- 
ing matter from distances <~ r m is responsible for bring- 
ing in positive angular momentum. At distances r ~ r m 
the angular momentum flux is carried mainly by the mat- 
ter. These results partially confirm the hypothesis of S94, 
OS95, Li et al. (1996) that matter should carry little an- 
gular momentum as it approaches the star. However, our 
results indicate that this does not lead to the "torqueless" 
accretion predicted by these authors, because the angular 
momentum flux carried by the matter is transferred to a 
flux carried by the magnetic field with decreasing distance 
from the star. The small twist of the magnetic field near 
the star, \B^\/B p <~ 0.1 (sec §6), is sufficient to carry the 
observed angular momentum transport flux (see discussion 
by Wang 1997). 

From our simulations we find that the magnetic field 
lines responsible for spin-up or spin-down of the star pass 
through the disk at distances r < r*. Consequently the 
spin evolution of the star depends on the location of the 
corotation radius r cor relative to r$. A star may spin- up, 
spin-down, or be in the regime of the "torqueless" accre- 
tion (see §4.2). For the case considered in this subsection 
of a slowly rotating star r cor = 3, the star spins-up. 

4. ACCRETION TO RAPIDLY ROTATING STARS 

We performed a set of simulations for different angular 
velocities of the star O*. Simulations were done for Type 
I initial conditions and with viscosity a = 0.01. The grid 
was Nr x Ng = 100 x 51 corresponding to the smaller 
region R max = 10. The fastest rotation considered was 
such that at r cor = 1, fi* = 1 which corresponds to a ro- 
tation period Po = 1.8 days for the T Tauri parameters of 
§2.2. (The conversion formulae for dimensional period is 
P* = 1.8/0* days). In this section we discuss the depen- 
dence of the spin-up/spin-down rate of the star on f2*. 

4.1. Rapidly Rotating Star: Spin-Down 

Here we discuss the case of a rapidly rotating star, 
f2* = 1, r cor = 1 (P* = 1.8 days), where r cor w r m . 
We observed that initially the accretion disk moved out- 
ward because the magnetic field of the star transferred its 
angular momentum to the disk. Later, for T > 3 the disk 
started to move inward towards the star, and finally a fun- 
nel flow formed (see Figure 11). In this case the thickness 
of the inner part of the disk is larger than that in the case 
of slowly rotating star, and matter accretes preferentially 
along the top layers of the disk. 

Figure 12 shows the observed fluxes of matter and angu- 
lar momenta to the star. Note, that the flux of an angular 
momentum is carried mainly by the magnetic field, Lf, 
which is negative so that the star spins down. Thus, we 



observe that as matter accretes to the star angular mo- 
mentum flows out from the star. As in the case of slowly 
rotating star, we observe a correlation between the matter 
and angular momentum fluxes. 

4.2. What is the Value of f2* for Torqueless Accretion? 

We arranged a set of simulations to test if a particu- 
lar value of O* gives "torqueless" accretion. According to 
GL79b, at some angular velocity of the star Cl cr it, posi- 
tive magnetic torque associated with the region r < r cor 
cancels the negative magnetic torque associated with the 
region r > r cor and the star can accrete without changing 
its angular momentum. 

We performed simulations at different f2* (correspond- 
ing to r cor = 1.2, 1.3, 1.4, 1.5, 1.7, 2, 3, and 10) and ob- 
served, that for fi* w 0.54 (r cor = 1.5) the angular momen- 
tum flux L wanders back and forth around zero (see Figure 
13). That is, it is positive for some time (<~ lOPo), then 
becomes negative for a similar length of time, and then be- 
come positive again. We conclude that this O* corresponds 
to torqueless accretion. The explanation for torqueless ac- 
cretion was proposed in the 1970s (e.g., GL79b). Namely, 
when more matter is supplied to the disk, the magneto- 
spheric radius r m moves closer to the star, and as a result 
there is a positive angular momentum flux to the star. 
When less matter accretes, the magnetospheric radius r m 
moves outward, and the flux of angular momentum is neg- 
ative. Figure 14 shows the evolution of the density and 
poloidal magnetic field for times corresponding to posi- 
tive angular momentum flux to the star (left panels), zero 
angular momentum flux (middle panels) and negative an- 
gular momentum flux (right panels). The variations of 
the accretion rate may be connected with accumulation of 
matter at r ~ rv where in this case 7V w 3 — 3.5 and 
subsequent accretion through instabilities. 

The ratio of the corotation radius r cor to the magneto- 
spheric radius r rn is of general interest. At the moments 
when the torque is zero, L = (see Figure 13), the mag- 
netospheric radius had the following values: at T = 24, 
r rn = 1; at T = 35, r m = 1.07; at T = 46, r m = 1.1. That 
is, for L = the ratio is: r cor /r m w 1.4 — 1.5. For sim- 
ilar simulations for faster rotation of the star f2» = 0.67 
(r CO r = 1-3), the angular momentum flux also wandered 
around zero value (see Figure 15), but most of time it 
was negative so that the star spins down. For simulations 
with more slowly rotating stars or larger corotation radii, 
r cor > 1.5, (r cor = 1.7, 2, 3, & 10), we observed that the 
star spins up. 

5. DEPENDENCE ON THE STELLAR MAGNETIC FIELD 



Fig. 15.— Same as on Figure 9, but f2* = 0.67, r cor = 1.3 (Po = 2.7 days for the T Tauri parameters of §2.2). The 
dot-dashed line corresponds to L m . The star spins down most of the time. 
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Fig. 16. — Disk-star interaction at different magnetic momenta of the star p after T = 10 rotations. Density (back- 
ground) changes from p = 2.5 (black) to p — 0.004 (white) Magnetic flux ^ (solid lines) is zero at panel a, and has limits: 
\I> = 0.1 — 0.6 on other panels. 



Fig. 17. — Close view of a quasi-stationary funnel flow calculated for Type I initial conditions. The color background 
represents the logarithm of the density. The density has a minimum value p — 0.004 (blue) and a maximum value p = 2 
(red). The lines are poloidal magnetic field lines or cquivalcntly ^(r,z) — const lines. The arrows represent the poloidal 
mass flux density pw. The bold dark line is the field line approximately in the middle of the funnel flow. The variation of 
different quantities along this field line are shown in subsequent figures. 



We did simulation runs for different magnetic moments, 
p = 0, 0.21, 0.42, 0.63 and 1.1, with other parameters 
fixed, fl* = 0.19 (r cor = 3), and a = 0.01. At T = the 
inner radius at of the disk was at r = 1 in all cases. The 
grid was similar to the one used in §4. 

Figure 16 shows result after T = 10 rotations. One can 
see that at smaller values of p, the inner radius of the disk 
and FF are settled closer to the star than in the case of 
larger p. In all cases, we observed that the base of the 
FF approximately coincides with the magnetosphere ra- 
dius r m , but at smaller p this radius is smaller in accord 
with equation (11). 

In case of even smaller p the direct accretion to the sur- 
face of the star is observed, like in simulations by HSM96 
and in some simulations by MS97. 

The cases with stronger magnetic moment p = 1.1, 0.63 
are more appropriate for explanation of T Tauri stars, be- 
cause the gap between the star and the disk is comparable 
or larger than the radius of the star. At smaller magnetic 
momenta p — 0.42, 0.21, the gap is too small. 

6. PHYSICS OF FUNNEL FLOWS 

Here, we discuss the physics of a typical funnel flow, 
the case of a slowly rotating star with r cor = 3 (§3) af- 
ter T = 30 rotations. The specific heat ratio is 7 = 5/3. 
Figure 17 shows an enlarged picture of the FF. We deter- 
mined the variation of different parameters as a function 
of distance (s) from the surface of the star along a fidu- 
cial poloidal field line through the middle of the FF. This 
field line is shown by the thick black line in Figure 17. 
Figures 18 - 21 show the variation of different parameters 
along this field line. Figure 18 shows that the density p 
and pressure P decrease along the FF, then increase again 
near the surface of the star. The temperature T initially 
varies slowly, but increases strongly approaching the star 
as a result of gas compression. 

Figure 19 shows the variation of different velocities along 
the FF. Matter flows along the FF with poloidal velocity 
v p . The velocity perpendicular to the poloidal magnetic 
field, v q , is very small, \v q \ « v p . The poloidal velocity v p 
increases and becomes larger than the slow magnetosonic 
velocity c sm at ~ 0.4 of the way. This point is expected 
to be closer to the disk in the case of cooler (and thinner) 
accretion disks (Koldoba et al. 2002). The flow become 
supersonic with the Mach number reaching M. s=s 3.6 at 
the surface of the star. Near the star, the poloidal velocity 



is close to the free-fall velocity: v p w Q.7vff. 

The poloidal velocity of the FF is much smaller than the 
Alfven velocity va- The ratio decreases from v p /va ~ 1/3 
in the middle of the FF, to v p /va ~ 0.07 at the star. 
Thus, the flow is strongly sub-Alfvenic. This validates the 
analysis of Koldoba et al. (2002) which proposed that the 
funnel flows were sub-Alfvenic. 

It is important to understand the forces which drive 
matter out of the equatorial plane of the disk into the 
funnel flow, and which forces accelerate the flow towards 
the star. The force per unit mass F along the magnetic 
field line is obtained by multiplying the Euler equation by 
Bp/ 1 Bp I . This gives 



to r sin x 



1 dp 

p ds 



9$ 
~ds~ 



1 d(rB^ 



8npr 2 



ds 



(12) 



(Ustyugova et al. 1999), where \ 1S the angle of the field 
line to the z— axis. The terms on the right-hand-side rep- 
resent the centrifugal force, the pressure gradient force, 
the gravitational force, and the magnetic force. 

Figure 20 shows forces in the outer two-thirds of the 
poloidal path of the FF. We see that the pressure gradient 
force is responsible for driving matter up and out of the 
disk and into the funnel flow, while the gravitational force 
is responsible for the acceleration of matter in the rest of 
the FF. Note that the magnetic force in equation (12) is 
small and acts in the opposite direction. Thus, we did not 
observe the "magnetic levitation" force predicted by Li & 
Wilson (1999). We observed a finite field above the 
disk in the FF (see Figure 21), but the magnetic force is 
much smaller than the other forces and it is not sufficient 
to levitate matter above the disk. 

Figure 21 shows that the angular velocity of the funnel 
flow decreases gradually from its value at the inner edge 
of the disk (r = Rd) to a small negative value at the star. 
The negative value of u) at the star is due to the fact that 
the magnetic field has a small inclination in the direction 
of rotation of the star. For this reason matter moving 
along the magnetic field line comes to the star with an 
angular velocity less than fi* (see also Ghosh et al. 1977; 
Muzerolle et al. 2001). 

Figure 21 shows the variation of the toroidal magnetic 
field along the funnel flow. The twisting of the dipolar field 
is relatively small with the "twist" 7^ = \B^\/B p < 0.15. 
At the star, 7^ w 0.04. This is in accord with estima- 
tions by Wang (1997), who has shown that the twist near 



Fig. 18.— Variation of the density p, pressure P and temperature T with distance s from the star's surface along the 
fiducial magnetic field line shown in Figure 17. 
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the star should be small, but not a zero. Away from the 
funnel flow \B^,\ is of the order of B p . We never observed 
conditions where >> B z (discussed, e.g., by GL79a 
and Agapitou & Papaloizou 2000). The reason is that for 
\B$ | >> B p , a torsional Alfven wave forms and propagates 
outward into the corona where it dissipates. 

In reality, matter may accrete through the magne- 
tosphere owing to 3D instabilities, for example, the 
Raylcigh- Taylor or Kelvin-Hclmholtz instability (Arons & 
Lea 1976a, b, Spruit & Taam 1990). Realistic investiga- 
tions of the role of these instabilities is possible with 3D 
MHD simulations. 

7. OPENING OF CORONAL MAGNETIC FIELD AND 
OUTFLOWS 

In our simulations we observed coronal magnetic field 
lines at large distances from the star opening due to the 
initial differential rotation. Closer to the star the field 
lines were mostly closed, contrary to the predictions by Aly 
(1986), LRBK95 (see also Newman, Newman, & Lovelace 
1992; Livio & Pringle 1992; Lynden-Bell & Boily 1994; 
Mikic & Linker 1994, Amari et al. 1996; FE99; FE00; 
Ustyugova et al 2000; Uzdenskyi, Konigl & Litwin 2002). 
The difference in behavior can be explained as follows. The 
analytical or semi-analytical models suppose that the in- 
ertia of the of the matter is very small, (3 = 8ttp/B 2 << 1, 
so that the plasma is in the "coronal" limit of Gold and 
Hoyle (1960). However, the inertia of the matter may be 
important. Romanova et al. (1998) simulated the dy- 
namics of magnetic loops in the corona of a differentially 
rotating disk where /3 ~ 1. They observed that in the mag- 
netically dominated regions (/3 < 1) the magnetic loops 
opened whereas in the matter-dominated regions (/3 > 1) 
the loops did not open. 

The energy-density of the dipole field decreases rapidly 
with distance (~ R~ 6 ) compared with the matter energy- 
density in our simulations which varies slowly within the 
computational region. For this reason the magnetically- 
dominated region is restricted to a relatively small volume, 
R < 1 — 1.5, around the star (see black line at Figure 5). 
To overcome this feature of the dipole field, GWB97 and 
GBW99 adopted a rapid fall-off of density in the corona, 
p ~ R~ 4 . For this dependence the region of small (3 is 
greatly expanded and opening of the magnetic field is fa- 
vored. 

In order to investigate possible stronger outflows, 
we chose parameters so as to increase the region of 
magnetically-dominated corona. The magnetic moment 
was increased 5 times, p = 5.3 (which corresponds to a 
magnetic field B = 5.3 kG at the surface of the star). The 
density in the corona was decreased 3 times, p c — 0.003. 
This led to a higher Alfven speed near the star, so that we 
increased the inner radius R m in = 0.8 to cover the region 
of high Alfven velocity, and increased the resolution of the 
grid to Nji x Ng = 120 x 71. We have f3 = l&tr = R d = l 
as in the other runs described in this work. The magnetic 
field is stronger and the initial density in the disk is about 



10 times higher than in other simulations. 

Figure 22 shows the results of the simulations. We ob- 
served that magnetic braking led to accretion of the inner 
part of the disk with velocity v r w 0.01 — 0.1. Thus matter 
in the disk concentrated in a dense ring with p ~ 10 — 14. 
The disk angular velocity is significantly smaller than the 
Keplerian value but larger than the angular velocity of 
the star, SI* = 0.19, which is a result of the interaction 
with the strong stellar magnetic field. A significant re- 
gion above the corona corotates with the disk. Analysis 
of forces, similar to that of §6 shows that in the region 
r < 2, the centrifugal force drives matter along the mag- 
netic field lines, but the force is not strong enough to open 
the magnetic field lines and matter accretes to the star 
forming multiple funnel flows (see Figure 22). However, 
the outer magnetic field lines, which are located at the 
edge of the magnetically-dominated region, become grad- 
ually open owing to combination of centrifugal and mag- 
netic forces. Matter accelerates along the region of open 
field lines up to velocities v ~ 0.5f eS c, where v esc is the 
escape velocity. We anticipate that if the density in the 
corona decreases faster than in our model, then matter will 
outflow from the region. The observed matter and angular 
momentum fluxes of outflowing matter are small. 

This test simulation provides evidence of how outflows 
can be formed from the disk in the quasi-stationary regime. 
Subsequent research is needed in this direction at different 
distributions of density of matter in the corona. 

8. APPLICATION TO T TAURI STARS 

Our simulations are well suited for modelling funnel 
flows to T Tauri stars because the entire flow can be well 
resolved. In the observed systems, the ratio of the inner 
radius of the disk where the FF begins to the stellar ra- 
dius is Rd/R* ~ 2 — 4. Such conditions are realized in our 
simulations. 

The detailed nature of the magnetic fields of T Tauri 
stars is not known, and the fields may differ significantly 
from a dipole field (e.g., Safier 1998). However, recent 
spectral studies of T Tauri stars show broad emission lines 
which give evidence of matter inflow with nearly free-fall 
velocities (e.g., Hartmann et al. 1994; Lamzin, et al. 1996; 
Muzerolle et al. 2001). Such velocities may be easily ex- 
plained in the model where matter inflows supersonically 
in a funnel flow along a dipolar field. 

In the following we perform numerical estimates for 
the funnel flow to a typical T Tauri star using conver- 
sion scheme of §2.1. We assume the inner radius of the 
disk is R = 2.8R*, the star's mass M . 8 = M/(0.8M Q ), 
its radius R 2 = R/(2R Q ), and its surface magnetic field 
£> 3 = £>/(10 3 G). Using these values we can derive refer- 
ence values for the different physical variables. For the flow 
velocity, v = {GM/R ) 1/2 w 1.63 x 10 7 M^R 2 1/2 cm/s. 
For the period of rotation at r = Rq, Pq = 2nR /v ~ 

1.78 M~j /2 R 3 2 /2 days. For the density, p = B%/v 2 w 
6.93 x 1Q- 12 B%R 2 Mq I g/cm~ 3 . For the accretion rate, 



Fig. 19. — Variation of different velocities with distance s along the fiducial magnetic field shown in Figure 17. v p is the 
poloidal velocity (bold line), v q is the velocity perpendicular to the magnetic field line, c sm and c/ m are the slow and fast 
magnetosonic velocities, Vff is a local free-fall velocity. The dashed line shows the Mach number M. 
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Fig. 20. — Forces parallel to the flow as a function of distance s from the star's surface along the fiducial field line in 
Figure 17. Only part (~ 2/3) of length of the flow is shown in order to resolve the forces acting near the disk. These 
forces are responsible for driving matter into the funnel flow. Fq is the gravitational force, Fq the centrifugal force, Fp 
the pressure gradient force, Fm the magnetic force, and Fr is the resultant force (bold line). 

Fig. 21. — Angular velocity u> and the "twist" \B^\jB v as a function of distance s from the star's surface along the bold 
field line in Figure 17. 

Fig. 22. — Density and poloidal magnetic field lines for a case with outflows. The magnetic moment of the star p = 5.3 
is 5 times larger than that used in most of our simulations. The density in the corona is 3 times smaller than in other 
simulations. Color background represents the density, white solid lines - contours of magnetic flux, black solid line is the 
line (3=1. Level lines of density are exponentially distributed between the minimum, p = 4x 10~ 4 , and the maximum, 
p = 17.6 values. The \& lines are exponentially distributed between "J" = 0.1 and "J" = 3. Arrows are velocities. The grid 
used was Nr x Ng = 110 x 71, and the viscosity parameter was a = 0.01. 



M = Po v Rl w 2.8 x \Q- 7 BlM~l /2 Rl /2 M /year. For 
the temperature of the corona, T w 1.07 x IO^Mq.sR^ 1 K. 

In the following we calculate the different physical vari- 
ables along the middle of the funnel flow shown in Fig- 
ure 17. From Figure 19 the matter reaches a maximum 
poloidal velocity Vp* w 277 km/s at the star's surface, 
which is 79% of the free-fall velocity Vff « 350 km/s. The 
funnel flow will of course give a distribution of velocities 
along the line of sight to a distant observer. The smallest 
line of sight velocities are for matter immediately above 
the disk where the speeds are ~ 10% of v p *. The line of 
sight velocity distribution is important for the interpreta- 
tion of observed spectral lines ( e.g., Folha and Emerson 
2001). The poloidal velocity is ~ 1/3 the Alfven velocity 
va above the disk (at s w 1) and decreases to v p rts va/15 
at the star's surface. The velocity across the magnetic 
field lines v q is very small (see Figure 19). Thus, matter 
flows along the magnetic field. The flow become supersonic 
above the disk and the Mach number reaches M. w 3.6 at 
the surface of the disk. 

The initial density of matter in the disk at Rq is p ~ 
6.9 x 10 -12 gcm~ 3 , and the density of matter in the FF 
(see Figure 18) varies in the range (1.2 — 9.2) x 10~ 12 cm~ 3 
which is typical for T Tauri stars. Our typical accretion 
rate to the star through the FF is M ~ (0.6 - 1.2) x 
10 -7 MQ/year, which agrees with estimates for T Tauri 
stars (Hartmann et al. 1998). The initial temperature in 
the disk is T w 9.6 x 10 3 K which is several times higher 
than expected in the disks of T Tauri stars at Rq. As men- 
tioned earlier this higher temperature is needed to give ad- 
equate resolution in the vertical direction through the disk. 
We observed that matter in the FF was heated by adia- 
batic compression to higher temperatures, T <~ 1.2 x 10 5 K. 

The fact that we use "free" boundary conditions at the 
surface of the star means that we cannot account for the 
stand-off shock at the stellar surface, surface heating, and 
other physical phenomena close to the star. Study of these 
processes can be done with a separate set of simulations. 

The investigated quasi-stationary funnel flows corre- 
spond to Classical T Tauri stars (CTTS) (see Hartman 
1998) where variability is not strong and possible out- 
flows are weak. Previous simulations of very fast accre- 
tion with subsequent strong outflows were done by HSM96, 
MS97, H97, GWB97, and GBW99, and these may corre- 



spond to earlier stages of evolution of young stellar ob- 
jects (YSOs of classes and I) which show evidence of 
episodic non-stationary accretion (e.g., FU Orionis type 
stars) and strong outflows (e.g. HH30). Non-stationary 
accretion may be initiated by thermal instability (Light- 
man & Eardlcy 1974) or by global magnetic instability of 
the disk (Lovelace, Romanova, & Newman 1994; Lovelace, 
Newman, & Romanova 1997), or by tidal interaction with 
a binary companion (Larson 2002). 

9. CONCLUSIONS 

9.1. Funnel Flows 

We have done a wide range of MHD simulations of disk 
accretion to a rotating aligned dipole in order to under- 
stand the different accretion phenomena. The simulations 
show that funnel flows (FF), where matter flows out of the 
disk plane and essentially free-falls along the stellar mag- 
netic field lines, are a robust feature of disk accretion to a 
dipole. Specifically, we find that 

(1) . The disk truncates and a funnel flow forms near 
the magnetosphere radius r m , where magnetic pressure of 
the dipole is comparable to the kinetic plus thermal pres- 
sure of the disk matter. 

(2) . The velocity of matter in the FF is much smaller 
than the Alfven velocity, |v| <~ (0.05 — 0.3)im, so that mat- 
ter flows along the magnetic field lines. The funnel flow 
accelerates and become supersonic. The Mach number is 
M. k> 3 — 4 at the surface of the star. At the star velocity 
is close to the free-fall velocity: v p w 0.7i>//. 

(3) . The angular velocity of the FF gradually varies 
from its value at the inner edge of the disk to the angular 
velocity of the star. The 'twist' of the magnetic field lines 
in the FF is small, \B^\lB v < 0.1, and it has a maximum 
approximately in the middle of the FF. 

(4) . The main forces which are responsible for drag- 
ging matter to the FF are matter pressure gradient force 
(near the disk) and gravitational force in the rest of the 
FF. Magnetic force is negligibly small. 

(5) . About 1/3 of the magnetic flux responsible for the 
spin- up/spin-down the star goes through the FF, while the 
remainder is above the FF. 

9.2. Disk- Star Interactions 
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Regarding the interaction between the disk and the star 
we find that 

(1) The magnetic field of the star influences the nearby 
regions of the disk inside a radius r br , while viscosity dom- 
inates at larger radii. The radius 7V depends on magnetic 
moment of the star /i and density in the disk. 

(2) . Inside the radius 7V the disk is strongly inhomo- 
geneous. The density is 2 — 3 times smaller than in the 
disk without magnetic field. Magnetically braked matter 
accumulates near magnetosphere and forms a dense ring 
and funnel flow. 

(3) . The star may spin- up or spin-down depending on 
the ratio of its rotation rate to the rotation rate at the 
inner radius of the disk. We find that "torqueless" ac- 
cretion is possible when r cor /r m w 1.5, where r cor is the 
corotation radius. 

(4) . At the star's surface, the angular momentum flux 
is transported mainly by the 'twist' of the magnetic field. 
Angular momentum carried by matter in the disk at <~ r m 
is transferred almost completely to the magnetic field at 
the star's surface. 

(5) . The coronal magnetic field is observed to open 
and close, but strong outflows were not observed for the 
considered parameters and quasi-equilibrium initial con- 
ditions. In the area of the disk where the field is strong, 



r < r$ « 0.5rj, r , the magnetic field lines tend to deceler- 
ate/accelerate rotation of the disk instead of being opened. 
Besides, opening of magnetic field loops is suppressed by 
matter dominated corona compared to GBW99 who ac- 
cepted dependence p ~ R~ 4 . 

(6). Strong outflows may be associated with strongly 
non-stationary accretion in the disk as observed in sim- 
ulations by HSM96, MS97, and H97 or when the disk is 
very dense and magnetic field lines are inclined to the disk 
as in simulations by FE99, FEOO. Alternatively the out- 
flows may occur in cases where the density of the corona 
decreases sufficiently rapidly with distance as in the sim- 
ulations by GWB97 and GBW99. Sporadic outflows can 
arise from a global magnetic instability of the disk (e.g., 
Lovelace et al. 1994, 1997, 2002), but they may be absent 
during the quiescent evolution of the disk-star system. 
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